Hyperfine Interactions in Graphene and Related Carbon Nanostructures 
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Hyperfine interactions, magnetic interactions between the spins of electrons and nuclei, in 
graphene and related carbon nanostructures are studied. By using a combination of accurate first 
principles calculations on graphene fragments and statistical analysis, I show that both isotropic 
and dipolar hyperfine interactions can be accurately described in terms of the local electron spin dis- 
tribution and atomic structure. A complete set of parameters describing the hyperfine interactions 
of 13 C and other nuclear spins at substitution impurities and edge terminations is determined. 
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Graphene and related carbon nanostructures are con- 
sidered as potential building blocks of future electronics, 
including spintronics [l[ and quantum information pro- 
cessing based on electron spins Q or nuclear spins 
Carbon nanostructures are attractive for these applica- 
tions because of the weak spin-orbit interaction in ma- 
terials made of light elements 0, Sj- Promising results 
for the spin-polarized current lifetimes in carbon nan- 
otubes [a 0, S| and graphene unambiguously confirm 
the potential of these materials. A number of quantum 
dot devices, components of solid-state quantum comput- 
ers, based on carbon nanostructures have been proposed 
recently 13, 11, I2> HI- Hyperfine interactions (HFLs), 
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the weak magnetic interactions between the spins of elec- 
trons and nuclei, become increasingly important on the 
nanoscale. In carbon nanostructures the interactions of 
electron spins with an ensemble of nuclear spins are ex- 
pected to be the leading contribution to the electron spin 
decoherence [!, 0, E3|. Minimizing HFIs is thus neces- 
sary for achieving longer electron spin coherence times 
[la ]. In some other instances the HFIs play an impor- 
tant role as a link between the spins of electrons and 
nuclei in the nanostructures 0, [H, IS EH underlying 
the implementations of quantum information processing 
involving nuclear spins. Probing HFIs with magnetic res- 
onance techniques also provides a wealth of information 
about structure and dynamics of carbon materials [lj|. 
A common understanding and an ability to control the 
HFIs are thus necessary for engineering future electronic 
devices based on graphene and related nanostructures. 

In this Letter, I study the hyperfine interactions in 
carbon nanostructures by using a combination of accu- 
rate first principles calculations on graphene fragments 
and statistical analysis. I show that the interaction of 
the conduction (low-energy) tt electron spins with nuclear 
spins can be described in terms of only the local (on-site 
and first-nearest-neighbor) ir electron spin distribution 
and the local atomic structure. The conduction electron 
spin distribution can be determined using simpler com- 
putational approaches (e.g. tight binding or analytical 
approximations [20L l2ll |) and tuned by tailoring nanos- 
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FIG. 1: (Color online). Projections of the spin-polarized con- 
duction electron density Pc( r ) ( a ) an d the total spin density 
p s (r) (b) on the plane of an electron-doped graphene frag- 
ment (in a.u. -2 ). The total spin density p s (r, z~0) (c) in the 
plane of nuclei (in a.u. -3 ) reflects the isotropic hyperfine field. 
Molecular framework is shown by black lines. 



tructure dimensions and applying external fields |22|, |23| | . 
The local nature of HFIs justifies the extension of my re- 
sults from small molecular models to extended systems. 
I further extend the considerations to curved topologies 
and to the presence of heteronuclei at impurities and 
boundaries. 

The all-electron density functional theory (DFT) cal- 
culations [24] ] were performed using a combination of the 
EPR-III Gaussian orbital basis set [25|| specially tailored 
for the calculations of hyperfine couplings and the B3LYP 
exchange-correlation hybrid density functional [26j . This 
computational protocol (see Ref. [27| for details) can be 
applied to molecules of limited size and predicts hyper- 
fine coupling constants (HFCCs) in excellent agreement 
with experimental results [27| . Spin-orbit and relativis- 
tic effects which are not important for the calculation of 
HFIs in light-element systems [28[ have been neglected. 
For a set of representative experimentally measured 13 C 
isotropic HFCCs of grapheme ion-radicals [13] our com- 
putations provide a mean absolute error of 1.1 MHz 
(«2% of the range of magnitudes), which justifies the 
use of calculated HFCCs as a reference. 

The effective spin-Hamiltonian of the HFI between the 
electron spin S and the nuclear spin I can be written as 



H=S A I, where the 3x3 HFI tensor A=A iso - 1 + T 
is usually decomposed into the scalar HFCC A iso and 
the traceless dipolar HFI tensor TT 30]. The HFI ten- 
sor A reflects the distribution of the electron spin den- 
sity p s (r)=pT(r)— pJ-(r) viewed from the position of the 
nucleus /. In carbon nanostructures the nuclear spins 
are those of the 13 C isotope natural abundance 

and can be atrificially changed; the dominant 12 C iso- 
tope has zero spin) and other elements originating from 
impurities and boundaries. The electron spin density 
p s (r) can be further decomposed into the contribution 
of half-populated conduction electron states lying close 
to the Fermi level (or singly occupied molecular orbitals 
in the molecular context), Pc( r )=J2 c IV' c ( r )| 2 ^0) an d the 
contribution of the fully populated valence states per- 
turbed by the exchange with spin-polarized conduction 
electrons, p s v (r)=J2v \ip v H r )\ 2 ~\' l P v H r )\ 2 ■ The crucial 
role of the exchange-polarization effect is illustrated with 
a model electron-doped hydrogen-terminated graphene 
fragment in the doublet spin state (Fig. [I}. While the 
projection of p s c {v) on the xy plane (Fig. QJi) is pos- 
itive everywhere and reveals an enhancement at the 
zig-zag edges, the projection of the total spin-density 
p s (r) (Fig. [TJb) is negative where p c (r) ^ s cl° se to zero. 
The isotropic (Fermi contact) HFCC is proportional 
to the total spin density at the position of nucleus /, 
Aiso={^Tt/'SS)l3 e (5Ngegip s {Yi) 1 where f3 e and (5n are the 
Bohr and nuclear magnetons, while g e and gi are the g- 
values of free electron and nucleus /, respectively. S is 
the maximum value of the electron spin projection. For 
the ideal graphene and planar sp 2 carbon nanostructures 
(all nuclei lie in the z~0 plane) p s c (z=0)=0 due to the p z 
symmetry of the conduction states. However, there is a 
contribution of the a symmetry valence states p s v (z=0)^0 
due to the exchange-polarization effect. For the model 
graphene fragment p s v (z=Q) (Fig. [T}:) shows an alternat- 
ing pattern with a relative dominance of the negative spin 
density. Since the a states are situated well above and 
well below the Fermi level in sp 2 carbon nanostructures, 
the valence exchange-polarization phenomenon exhibits 
the property of locality. This property was exploited by 
Karplus and Fraenkel almost 50 years ago to describe 
the isotropic 13 C HFCCs in conjugated organic radicals 
[3lj ]. The main contribution to the hyperfine anisotropy 
originates from the total spin population n of the on-site 
p z atomic orbital, which also incorporates the contribu- 
tion of exchange-polarized valence states. Assuming a 
local axial symmetry, T can be written as a diagonal 
matrix with elements T zz /2=— T xx =— T yy =A ( n p , where 
Adip=(l/5S)(3 e (3Ng e gin(l/r2 P ) {t2 P is the distance of the 
carbon 2p electron to nucleus). 

The HFIs were calculated for a set of 12 (~1 nm 
size) electron- and hole-doped planar hydrogen termi- 
nated graphene fragments (Fig. [5^) in the spin-doublet 
ground states. This provides overall statistics for 206 in- 
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FIG. 2: (Color online) (a) Set of graphene fragments used 
in the present calculations. Fitted Fermi contact (A iso ) and 
dipolar (Adi P ) 13 C HFCCs vs the corresponding values, Ai SO 
(b) and Adi P (c), calculated from first principles. The values 
for inner (3 carbon NNs) and boundary (2 carbon NNs) atoms 
are shown as red and blue dots, respectively. 



TABLE I: Parameters (in MHz) fitted to the results of calcu- 
lations of a set of nanographite molecules. 
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equivalent 13 C HFCCs. The calculated Ai SO and A<n p 
values are fitted to the extended form of the Karplus- 
Fraenkel expression 

A = aj {l+ J2 bjArA^ + c (1 + dAr^nf, (1) 
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where the two terms account for the contributions of the 
on-site and nearest neighbor (NN) conduction electron 
spin populations per unpaired electron, n c and n?, re- 
spectively, calculated from first principles. The on-site 
coefficients a,j and bj are distinguished for the cases of C 
atoms with 3 carbon NNs (j=3) and the boundary atoms 
with 2 carbon NNs (j=2). The C-C bond length effects 
are encountered through the coefficients bj and d with 
^■ r i =r i~ r o being the deviation of the bond length from 
the value for the ideal graphene, ro=1.42 A. Only statisti- 
cally significant local properties were included in the lin- 
ear expression ((T|) . The results of the regressions are sum- 
marized in Tab.[fl(l MHz=4.136x 10~ 3 (usY). Fig. IU[b,c) 
shows the fitted (using expr. ([1]) and regression parame- 
ters) values Ai SO (Adi P ) versus the calculated Ai SO (Adi P ) 
values. Regressions to the linear expression (QJ pro- 
vide accurate estimations (root-mean-square-errors are 
1.7 MHz and 1.2 MHz for Ai SO and Adi P , respectively). 
The calculated isotropic HFCCs span about the same 
range of magnitudes (—16.7 MHz<^4i SO <21.1 MHz) as 
the dipolar HFCCs (-5.9 MHz<A dlp <22.7 MHz). The 
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HFCCs of boundary atoms tend to be larger due to the 
fact that low-energy states localize at the zigzag graphene 
edges [2(| . The on-site and the NN exchange-polarization 
effects have competitive character (et 3 /c«— 2) in the case 
of isotropic HFCCs. Our calculations predict «50% 
larger values for the parameters a2, 03 and c for Ai so 
compared to those obtained by Karplus and Fracnkcl 
in their early studies of HFCCs in molecular radicals 
(a 2 =99.8 MHz, a 3 =85.5 MHz and c=-39 MHz) 
This difference can be explained by the incorporation 
(via DFT) of the electron correlation effects in our calcu- 
lations and to the local atomic structure of the graphene 
lattice. Both Ai SO and A<n p show a tendency to enhance 
the on-site and to weaken the NN contributions with 
the increase of C-C bond lengths. The dipolar HFCC 
is mostly influenced by the on-site contribution of the 
half-populated conduction state and the NN exchange- 
polarization effect is weaker in this case (03 /c~— 7). 
When compared to typical solid state environments based 
on heavier elements, the 13 C HFCCs in graphene and re- 
lated nanostructures are weaker (e.g. 117 MHz 31 P Fermi 
contact HFCC for the P shallow donor in Si 32j) and 
more anisotropic. 

The graphene honeycomb lattice is a bipartite lat- 
tice, i.e. it can be partitioned into two complementary 
sublattices A and B. I discuss the HFIs for the three 
general cases of conduction electron spin distributions 
over the sublattices: (i) ferromagnetic n A =n c B >0; (ii) 
fcrrimagnetic n A >0 and n c B —Q and (hi) antiferromag- 
netic n A =—n B >0 (see Tab. HTJ. The first case can be 
physically realized upon the uniform magnetization of 
the system with equivalent A and B sublattices, e.g. 
by applying an external magnetic field. The negative 
Ai SO ——AA MHz is small due to the partial compensation 
of the on-site and the NN exchange-polarization effects. 
This value is consistent with the values derived from the 
experimental 13 C Knight shifts in graphite intercalates 
(-25 MHz<A J;5O <-50 MHz) ^ and with the calcu- 
lated isotropic Knight shifts in metallic carbon nanotubes 
34j . The ferrimagnetic case with the conduction state 
distributed over the atoms of only one sublattice (A) is 
physically realized at the zigzag edges [20| and around 
single-atom point defects in sublattice B [35(. Consid- 
erable alternating Fermi contact and dipolar HFCCs are 
predicted in this case. An antiferromagnetic pattern can 
be realized in the case of heavily disordered systems with 
localized defect and edge states in both sublattices [36| . 
The magnitudes of HFIs are minimized and maximized in 
the cases of ferromagnetic and antiferromagnetic electron 
spin distributions, respectively. 

Many carbon nanostructures of reduced dimension- 
ality (e.g. nanotubes and fullerenes) represent non- 
planar topologies. Local curvatures lead to the sp 2 ~sp 3 
rehybridization of carbon atoms and enable a Fermi 
contact interaction involving the low-energy 7r electron 
spins 19]. This results in a positive contribution of 



TABLE II: Hyperfine coupling constants for three general 
cases of spin populations n c of the carbon atoms of A and 
B sublattices of graphene. 





A iso {A)/n c 


A tso (B)/n c 
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n c A >0; n c B =0 


128 
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131 


-58 


n c A =~n c B >0 


300 


-300 


189 
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TABLE III: Parameters (in MHz) describing the HFIs of nu- 
clei of susbtitutional impurities ( n B and 14 N), monoatomic 
functional groups ( X H and 19 F) and rehybridized (sp 3 ) carbon 
atoms ( 13 C) at the edges. 
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Adip 

a c 
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subst. impurity 


43 


-31 


60 6 


14 N 


subst. impurity 


150 


-22 


130 -10 


l R 


C sp 2 edge 


-119 


22 




19 F 


C 3p 2 edge 


240 


-40 




I H 


C 3p 3 edge 




350 




19 p 


C 3p a edge 




750 




13 c 


C 3p3 edge 




-68 





the 7r states unless n c is close to zero: a contribu- 
tion due to the NN exchange-polarization effect is neg- 
ative in this case. The degree of rehybridization 777 
of the 7T states (s m p ) c an be described using a local 
bond angles analysis [37[. For the case of large curva- 
ture radii the original expression for m can be refor- 
mulated in a more convenient form, m=d 2 , c j '8(1 /Ri + 
I/R2) 2 , where d cc is the C-C distance, R\ and R2 are 
the principal curvature radii. The curvature-induced 
contribution to the Fermi contact 13 C HFCC is then 
A%™={4Tr/3S)p e f3 N g e9l nmc!)l(0), where 2s (O) is the 
magnitude of the carbon atomic 2s wavefunction at the 
point of nucleus ((87r/3)/3 e /3jvg e 3/^ s (0)«3.5xl0 3 MHz). 
The curvature-induced direct coupling becomes signifi- 
cant (m>10 -2 ) only in ultranarrow carbon nanotubes 
(d<l nm) and fullerenes. 

Since the natural abundance of the "HFI-active" 13 C 
isotope is small («1%), consideration of the nuclei of 
other elements is important for a complete description 
of HFIs in carbon nanostructures. The common substi- 
tution impurities are boron and nitrogen with all natu- 
ral isotopes having nuclear spins. Graphene edges can 
be terminated by hydrogen and fluorine atoms with 
both *H and 19 F spin-1/2 nuclei (99.9885% and 100% 
natural abundance, respectively) having high (/-values 
(,g( 1 H)/.g( 13 C)« ff ( 19 F)/ 5 ( 13 C)«4). I consider HFIs in a 
reduced set of molecular fragments (only 3- and 4-ring 
structures included) with impurities and edge functional- 
izations in all possible positions. The calculated HFCCs 
have been fitted to the Karplus- Fraenkel relation, with no 
Ar terms included (Tab. HTl]) . The variations of the local 
charge density of states in the vicinity of impurities does 
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not have any significant influence on HFIs. Both Fermi 
contact and dipolar HFCCs of the impurity nuclear spins 
show a monotonic increase along the 11 B- 13 C- 14 N series 
when compared to the results for 13 C HFCCs (Tab. H|. 
The NN relative exchange-polarization effects (a/c ratio) 
on the Fermi contacts HFCCs tend to decrease along the 
series. While the HFIs of the nuclear spins in substitution 
impurities are highly anisotropic, the hypcrfinc couplings 
of the edge nuclei show small anisotropy due to the sp 3 
character of bonding. When : H and 19 F edge nuclei are 
bound to the C sp 2 atoms, the isotropic HFCCs are of 
the same order of magnitude as those of the 13 C spins 
in the graphene lattice. The influence of the NN carbon 
atoms (second NNs to the terminating atom) is very sim- 
ilar for 1 H and 19 F nuclei and smaller than in the case 
of 13 C HFCCs (a H /c H ^a F /c F w~6). The spin polar- 
ization effect on 19 F HFCCs is stronger and of opposite 
sign compared to that of protons (of / <ih~cf / ch~— 2). 
When edge atoms are bound to the rehybridized (sp 3 ) 
carbon atoms, n c is zero but the NN contribution is sig- 
nificantly enhanced. The NN contribution to the 13 C 
hyperfine coupling of the sp 3 edge carbon atom itself 
(c=— 68 MHz) has a similar magnitude as that of the 
sp 2 edge atoms (c=— 57 MHz). HFIs with the boundary 
spins (H-terminated edges are often obtained in experi- 
ments [381 ]) have to be taken into account when designing 
carbon-based nanoscale devices for spintronics or quan- 
tum computing. A chemical modification of the graphene 
edges (e.g. substitution of the hydrogen atoms by alkyl- 
groups) can be suggested to reduce electron spin deco- 
herence effects from the HFIs with boundary spins. 

In conclusion, the results of first principles calculations 
show that the hyperfine interactions in graphene and re- 
lated nanostructures are defined by the local distribution 
of the conduction electron spins and by the local atomic 
structure. A complete set of parameters describing the 
hyperfine interactions was determined for the 13 C and 
other common nuclear spins. These results will permit 
control of the magnetic interactions between the spins of 
electrons and nuclei by tailoring the chemical and isotopic 
compositions, local atomic structures, and strain fields 
in sp 2 carbon nanostructures. Some practical recipes for 
minimizing interactions with nuclear spins are given. 
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